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In this work we investigate and parameterize the amount and angular distri- 
bution of Cherenkov photons which are generated by electro-magnetic cascades 
in water or ice. We simulate electromagnetic cascades with Geant-4 for primary 
electrons, positrons and photons with energies ranging from 1 GeV to 10 TeV. We 
parameterize the total Cherenkov light yield as a function of energy, the longitu- 
dinal evolution of the Cherenkov emission along the cascade-axis and the angular 
distribution of photons. Furthermore, we investigate the fluctuations of the total 
light yield, the fluctuations in azimuth and changes of the emission with increasing 
age of the cascade. 
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1. Introduction 

High energy neutrino telescopes such as IceCube, Baikal or Antares HI ID 13 
detect Cherenkov light from charged particles in natural media like water or ice. 
Cherenkov light is produced when these particles propagate through the medium 
with a speed faster than the phase velocity of light v > Cmed - c/n. In ice and water 
the refraction index n is typically « « 1.33 HIH. Hence, the Cherenkov threshold 
is given by /3 = ^ which corresponds to a minimum kinetic energy of 
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For electrons this is Ec ~ 0.26 MeV in water and ice. The number of emitted pho- 
tons per unit track x and wavelength interval is given by the Frank-Tamm formula 

mm 

(fN Inaz^ J 

(2) y-TT = ■ sin\ec) 

clxclA 

Here 9c is the Cherenkov angle. This is the opening angle of a cone into which the 
photons are emitted 

(3) cos(0,)-^ 

np 

A relativistic track (yS = 1) in water or ice ~ 1.33) produces about A/^o ~ 250cm~^ 
optical photons in a wavelength interval between 300 nm and 500 nm, which is a 
typical sensitive region of photo-detectors, e.g. |j8l, used in the aforementioned 
neutrino telescopes. The Cherenkov angle for a relativistic track (yS = 1) in ice is 
9c ft - arccos(l/n) w 41°. 

A large fraction of detected Cherenkov photons in high energy neutrino tele- 
scopes originates from electromagnetic cascades. These are initiated by a high 
energy electromagnetic particle which produces secondary particles by subsequent 
bremsstrahlung and pair-production processes ||6l [3 . The primary particle can 
originate from radiative energy losses of a high energy muon (bremsstrahlung and 
pair-production), or from the decay of 2y in hadronic cascades, or be a high 
energy electron from a charged current interaction of an electron neutrino. 



Figure 1: A simulated electromagnetic cascade. A primary electron of 100 GeV 
has been injected at the left pointing towards the right. Shown are all generated 
charged secondary particles (red for negative and blue for positive charge) as the 
result of a Geant4 simulation. Neutral particles, like photons are not shown. 
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An example of a simulated cascade is shown in figure [T] Each particle in 
the cascade produces Cherenkov light according to eqj2] and eqj3} if its energy 
is above the Cherenkov threshold eq{T] Due to multiple interactions and scattering 
the directions of the particles in the cascade differ from that of the original primary 
particle and a broad angular distribution of emitted Cherenkov photons is expected. 

The characteristic scale length for the development of an electromagnetic cas- 
cade is given by the radiation length Xq |6]. It is about Xajce * 39.75 cm and 



Xo,watei- * 36.08 cm as determined by Geant4 for the configuration listed in Ap 



pendix A The length along the shower axis z is usually expressed by the dimen- 
sionless shower depth 

(4) t = z/Xo . 

The length of the cascade increases typically logarithmically with the ratio of pri- 
mary energy Eq and critical energy Ecnt- The critical energy is the energy above 
which radiative process dominate the energy loss of electrons. The values for ice 
and water obtained from IM are E""-., . = 78.60 MeV, . = 76.51 MeV and 

" " cnt,ice ' cnt,ice 

, = 78.33 MeV, . ^76.24 MeV. 

crit, water ^ -^-^ > water,ice 

The physical length of a shower is typically less than 10 m. This is short, 
compared to the scale of neutrino telescopes and the full Cherenkov light is created 
locally and expands with time as an almost spherical shell with a characteristic 
angular distribution of the intensity. 

Due to the large number of particles the full tracking of each particle in Monte- 
Carlo simulations of cascades in neutrino telescopes is very time consuming. How- 
ever, the development of electromagnetic cascades is very regular because fluctu- 
ations are statistically suppressed by the large number of interactions and large 
number of involved particles. Hence, it can be well approximated by the average 
development. Therefore, for the simulation of data in neutrino telescopes, the aver- 
age Cherenkov light output can be parameterized, e.g. as done in |[m[T2l[T3l[T4ll . 

This work follows up the work in ffTTl which was based on Geant3. 16 |[T5l with 
a more precise calculation of the total Cherenkov-light yield and its angular distri- 
bution based on Geant4 [.16j . For different primary energies and primary particles, 
we investigate the velocity distribution and the directional distribution of particles 
and the longitudinal development of the cascade. We present a parameterization 
of the Cherenkov light yield and investigate its fluctuations as well as variations 
of the azimuthal symmetry of the cascade. We also present a parameterization of 
the angular distribution of Cherenkov photons and investigate variations of this 
distribution during the development of the cascade. The results are compared to 

uniiniinKii. 
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2. Simulation method 



The calculation of this work follows largely the strategy described in ifTTIl . 
We use the Geant4 (GEometry ANd Tracking) toolkit to track the particles in the 
cascade through the medium ice or water lll6l . The used media properties are 
given in Appendix A| Unless noted otherwise, we used an index of refraction of 
n = 1.33 and a density of p,ce = 0.91 g/cm^. Note, that these values slightly de- 



viate from the values in lUOl : p,ce = 0.918 g/cm^ and nice = 1.31 g/cm"* and the 
value Pice = 0.9216 g/cm^ at the center of IceCube ifTSl . This introduces a small 
systematic uncertainty of about 1%, which can be corrected for by rescaling our 
results to the correct density. 



— primary 

— charged particle 
■■- neutral particle 



60 m 



40m 





Figure 2: Geometry of the simulation and method of the calculation 

The simulation principle of this work is illustrated in figure |2] The medium is 
contained in a cylindrical volume of 30 m radius and 40 m height. The dimensions 
are chosen such that all secondary particles are well confined within the geometry 
and fully tracked. The primary particle e- or y is injected at the bottom center into 
this volume with its initial momentum pointing into positive z-direction. The par- 
ticles are propagated through the simulation and secondary particles are created, 
which again can produce further particles. Each step between two interactions 
corresponds to a track segment for which the energy and direction are assumed 
constant. For each track segment / we store the length the Lorentz factor the 
z-position z.i and the direction a, with respect to the z-axis. The azimuth angle O 



discussed in section 3.6 corresponds to the rotation angle in the x-y plane. Sum- 



ming over all track segments allows to calculate the Cherenkov photon yield and 
the corresponding angular distribution. 

For the here described simulation it is important to simulate all particles with 
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energies above the Cherenkov threshold. Details on the simulated physics pro- 



cesses are given in Appendix A In Geant4 some electromagnetic processes re- 
quire production thresholds to avoid infrared divergences fT9l. These production 
thresholds are specified as a cut-in-range threshold, using the SetCutsC) method 
of G4VUserPhysicsList. Here, particles are tracked if their mean expected range 
is larger than this cut-in-range threshold. For each material and particle type, this 
cut-in-range is transformed into a corresponding energy threshold. Here, a cut- 
in-range of lOOyum is chosen. This corresponds to a kinetic energy threshold 
of Ecut,e* ~ 80keV for electrons, which is well below the Cherenkov threshold 
Ec^e^ ~ 264 keV. Once produced, all secondary particles are tracked until they 
stop. In order to increase the computing performance, a single scattering process 
of a particle does not correspond to an individual track segment but multiple scat- 
tering processes are simulated as one step. 

We perform simulations up to a maximum energy of the primary particle of 
10 TeV to which electromagnetic cross sections in GEANT-4 are valid. Our results 
could be extrapolated beyond this limit, however, at high energies an additional 
effect, the LPM effect, is expected to set in. This effect describes significantly re- 
duced electromagnetic cross sections and the longitudinal development of such cas- 
cades would become strongly elongated such that our parameterization approach 
is not valid. 

For p = I, the number of emitted Cherenkov photons is proportional to the 
length of the track and can be calculated using equation |2] For yS < 1 the photon 
yield is smaller and proportional to the factor 

1 



(5) sin^iOc) = 1 - cos^(6/c) = 1 - 
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In order to properly account for this smaller yield, the length of each track segment 
/ is scaled with the Frank- Tamm factor 

sin^C^c) 9 1 

(6) 1= . ' with sm\ecfi) = l-^ 

The value / thus corresponds to the equivalent length of a relativistic track with the 
same photon yield as the track length /. The use of the equivalent length / instead of 
an explicit calculation of photons has the advantage that the here presented results 
can be rescaled to slightly different indices of refraction, and are independent on 
the assumed wavelength interval of the considered photo-detector. 

The angular distribution of the Cherenkov photons is calculated with the method 
introduced in [17 1. In this method the distribution of track length as a function of 
the directional angle a with respect to the z-axis (zenith) and the velocity yS can 
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be transformed to a zenith distribution of emitted Cherenkov photons. The prereq- 
uisite for the applicability of that method is a high statistics of tracks which are 
distributed uniformly in azimuth. 



3. Results 

3.1. Velocity distribution of shower particles 

Figure[3]shows the distribution of track length in the cascade versus the Lorentz 
factor yS for different primary energies Eq. The shape is remarkable constant for 
different Eq while the total normalization is proportional to Eq. The distribution 
also does not depend on the type of the primary particle. The difference between 
the physical track-length I and the effective track length / becomes particularly 
obvious close to the Cherenkov threshold ji w 0.752. Unless noted otherwise, we 
will use in the following the effective track length / instead of the physical track 
length I. 

3.2. Total light yield 

Figure |4] (top) shows as an example the distribution of the total summed ef- 
fective track length for repeated simulations of a primary electron of 1 TeV. The 
distribution can be well described by a Gaussian, which is fit to the data. 

The mean expectation and the standard deviation from Gaussian fits to distri- 
butions for different primary particles and different primary energy Eq are shown 
in figure [4] (middle) and (bottom). The data is fit with a power-law 

(7) t{Eo) = a-El , o-,<£'o) = a ■ E^^ 

In all fits the parameter p is found to be consistent with 1 at the level 10"^ indicating 
a very good linear relation between the total / and £"9. Also the coefficients a agree 
within 10~^ for different primary particles. The detailed results of these fits are 



given in table|B.3|in Appendix B As a result we obtain an energy scale parameter 



which relates linearly the total Cherenkov light yield with the primary energy 

(8) a w 532.1 ± 10~^ cmGeV"^ 

The observed value a = 532.1 cm GeV^' is slightly larger than the value a = 
521 cm GeV~' in d which was also obtained in Geant4 simulations of ice. The 
origin of this 2% difference is not obvious. It could be either related to difference 
in the versions of Geant4 but also to different configurations which is not given in 
|[T2| . e.g. a slight difference of the assumed index of refraction. 
If rescaled to the density of water by the relation 

(9) Pwater ' (^water ~ Pice ' (^ice 
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Figure 3: Velocity distribution of track length. Shown is the differential distribution 
of summed track length per shower versus the Lorentz factor ji for bins of 0.002 
in p. The top figures show the distributions of physical length / for the shower 
from primary positrons of different energies. The middle figure shows the same 
distribution for the track length / which has been weighted with the Frank-Tamm 
factor, eqj6] The bottom figure shows the distributions of / for different primary 
particles: e'^ ,e~ ,y for the primary energy £"0 = 1 TeV. 
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Figure 4: Total amount of Cherenkov-light-radiating track length I, including the 
Frank-Tamm factor. Top: distribution of / for 5 • 10^ simulated primary electrons of 
Eq = I TeV. A Gaussian-distribution is fit to the data. Middle: /(^o) as a function 
of the primary energy Eq resulting from Gaussian fits. A power-law function eqjV] 
is fit to the data. Bottom: Standard deviation o-f(Eo) Resulting from Gaussian fits. 
A power-law function like eqlTlis fit to the data. 
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we obtain the value abater * 484 cm GeV-i. This is significantly larger than the 
values 437cmGeV-i found in [11] for n = 1.33 and 466cm GeV"^ in Oil for 
n = 1.35. 

As expected for an increasing number of particles, the size of fluctuations of 
the total track length cTf increases oc ^/Eo and the value /? = 0.5 is fixed for the 
fit. Hence, the relative size of fluctuation decrease oc 1 / ^/Eo with higher primary 
energy Eq (see figure |4]). With the values in table B.3| in Appendix B the relation 



(10) 

is found. 



(Tf 



0.0108 



IGeV 



Eo 



3.3. Longitudinal cascade development 




distance z ccm] distance z [crrt] 



Figure 5: Longitudinal shower profiles as a function of the length z along the 
shower axis. Shown is the track length distribution relative to the total length Zq 
of the cascade. The left figure shows the result for initial y and different primary 
energies Eq. The right figure shows the result for different primary particles e-, y 
and a primary energy Eq = 100 GeV. All track segments have been weighted with 
the Frank-Tamm factor. 



The longitudinal profiles of the track length dl/dz along the axis of the cascade 
are shown in figure |5] As expected from a simple Heitler model BO, the depth 
of the shower maximum Zmax scales logarithmically with the primary energy. The 
distributions are almost identical for and e~. However, for an initial photon the 
depth of the shower maximum is about one radiation length deeper. 

The longitudinal shower profile can be parameterized with a gamma distribu- 
tion ^ 

dz (btr^ 

(11) /tot • — = b , 
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distance z [cm] 



Figure 6: Example of the parameterization of the longitudinal shower profile with 
eq 1 1 for a positron with a primary energy Eq = lOTeV. 



Here, t is the shower depth t = z/Xq, a and b are characteristic dimensionless 
constants of the cascade [6]. An example fit is shown in figure [6] The results of all 
fits for difi'erent Eq are given in table |B .71 in [Appendix B 




energy [GeV 



Figure 7: Fit parameters of the longitudinal profile a, b versus initial energy. 



The energy dependence of the fit parameters a and b is shown in figure [7] is 
found constant and does not depend on the particle type, while the parameter a can 
be described with an logarithmic increase 



(12) 



« = -+^-logio(Yfv) 



It is slightly larger for y than for e-. The parameterization results for b and a, fi are 
given in table B.4 in [Appendix B 
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Figure 8: Maximum of the longitudinal Cherenkov-radiating track-length profile 
as a function of the initial energy. The markers represent the calculated values for 
tmax and solid lines are for visual guidance. The dotted/dashed lines show equation 
([T4]) with the parameters from ||6l (PDG)/ |[20ll (Grindhammer) respectively. 



Figure [8] shows the shower maximum of the Cherenkov-radiating track length 
as a function of the initial energy. It has been calculated from the longitudinal 
shower profiles with the formula 



(13) 



Additionally shown is the maximum of the longitudinal energy deposition based 
on the parameterization in ||6l 



(14) 



tmax = Iny + Cj, i = e,y, 



with y = Eo/Ecrit- In [6] for electron- or positron-induced cascades the values 
Cg'-"^ = -0.5 and for photon-induced cascades C™*^ = -1-0.5 are given. These 
values are based on simulations with EGS4 up to an energy of 100 GeV for nuclei 
heavier than carbon. Up to that energy the slope of our result agrees but our values 
are offset by about -0.5. Above 100 GeV we also deviate in slope. 
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In contrary, EOl gives the value C^"^ = -0.858 for electron- or positron- 
induced cascades. This parameter was also obtained from fits to the longitudinal 
energy deposition profiles for elements ranging from carbon to uranium at ener- 
gies from 1 GeV to 100 GeV. The simulations were performed with GeantS. The 
value Cy^' is not explicitly stated. Assuming that the difference of the maxima 
of photon- and electron-induced cascades is about one radiation length we obtain 
C^'- « C2'- + 1 - +0.142. 

Here, the simulations are performed with Geant4 and ice is used as the detector 
material. It can be seen that our results agree much better with [20|. Nevertheless, 
deviations appear for larger energies, where the here obtained energy dependence 
increases less than logarithmically. This effect is stronger for photon-induced cas- 
cades than for electron-induced cascades. 



3.4. Angular distribution of tracks 




20 40 60 80 100 120 140 160 180 20 40 60 80 100 120 140 160 180 

zenith angle a [deg] zenith angle ct [deg] 



Figure 9: Distribution of the effective track length / per shower as a function of 
the inclination angle a of the primary particle's direction. The left figure shows 
the result for a primary positron and different Eq and the right figure the results for 
£0 = 1 TeV and different primary particles. 



Relevant for the angular distributions of Cherenkov photons is the angular dis- 
tribution of secondary tracks in the cascade and their velocity . Here, a is the 
polar angle of the track with respect to the z axis and the Lorentz-factor. 
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The distribution ^ is shown in figure 



and the distribution -f-jn in figure 



dad/3 ^ 

Most particles are produced in forwardHirection with a velocity j6 close to c. 
can be seen that the normalization of the o'-distribution changes with energy but the 
shape does not. The angular distribution does not change with for different primary 
particles. 
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velocity p 

Figure 10: Density distribution of the effective track length / versus the inclina- 
tion angle a and Lorentz factor /3 for a 1 TeV shower. The vertical color codes 
corresponds to the histogrammed lengths I per shower. 
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Figure 1 1 : The angular distribution of secondary Cherenkov photons for different 
Eq (left) and different primary particles (right). Shown are the normalized angular 
distribution per photon and steradian. 



3.5. Angular distribution of Cherenkov light 

The angular distribution of Cherenkov photons is calculated with the method 



described in lfT7l . Figure Tlj shows example distributions ^ versus the polar angle 
(p with the z-axi^ The distribution corresponds to the track length that produces 
an equivalent total Cherenkov light yield. 

A broad distribution with a clearly pronounced Cherenkov peak is visible. As 



'Note that the definition of is identical to the previously defined angle a. However we use a 
difl^erent symbol to indicate the difi'erence of photons and tracks. 
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expected from the results in section 3.4 the shape of the distribution is unchanged 
for different primary energies and different primary particles. 




Figure 12: Example of the parameterization of the angular distribution for a Eq = 
100 GeV positron. The parameters PltoP4 correspond to the parameters a to din 
eqfT^ 



The angular distributions are parameterized with a simple function 



(15) 



ae 



A typical fit is shown in figure 12 



The fit parameters for different energies are given in table B.9 in Appendix B 



They are found to be very similar and constant with energy. We conclude that the 
angular distribution of Cherenkov photons can be described with the above formula 



and the averaged parameters given in table B.5 in Appendix B 



3.6. Fluctuations in azimuth 
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Figure 13: Example of the azimuthal distribution of four individual showers of 
electrons with a primary energies £"0 = 10 GeV (left) and £"0 - 1 TeV (right). The 
relative effective track length per 30° -interval is shown. 
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An important prerequisite for the here used calculation of the angular distri- 
bution of Cherenkov photons is the assumed symmetry in azimuthal angle O of 
the distribution of track direction^ in the plane around the shower axis [17|. As 
an example, figure 13 shows azimuthal distributions of / for four individual show- 
ers, each for two different energies. Differences originate from fluctuations in the 
shower development. Correspondingly, the relative size of fluctuations strongly 
decreases for larger particle numbers in higher energy showers and the distribution 
becomes almost flat in azimuth. 

The effect of this asymmetry can be quantified by aligning all simulated show- 
ers in azimuth in order to account for their random orientation. To obtain an av- 
eraged azimuthal distribution the bin contents for each shower are added with the 
maximum bin aligned and the azimuthal orientation is defined according to the di- 



rection of the second highest bin. The results are shown in figure 14 The mean 
total amplitude of angular fluctuations in azimuth can be as large as ±11% for 
10 GeV but decreases approximately with the square root of the primary energy to 
less than ± 1 . 1 % for £0 = 1 TeV. However, as indicated by the error bars, the indi- 
vidual bin fluctuations are of the same order of magnitude as the mean amplitude 
of the asymmetry. 

The amplitude of the asymmetry is fit with the paraboUc function 



(16) 
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l+bUcp-cPoY 



4>of + 0o) 



The parameter b describes the vertical compression and (po the position of the min- 
imum. The third term accounts for the normalization. The angle is used in units 
of radians and the bin size of the histogram is AO = 30° = 0.523 rad. The results 
of all fits are given in table |B.10| in [Appendix B 



The parameter 0o is found roughly constant. It differs slightly from n because 
the preferential direction of the second largest bin leads to an angular bias. The 
energy dependency is shown in figure [14] (bottom right). The amplitude coefficient 
b is fit with 

(17) b^-^ 



as expected from the correspondingly increased number of particles in the cascade. 



The results of these fits are summarized in table B.l 1 in appendix Appendix B 



'Note, that the azimuthal angle <t> is diiferent from the earlier defined angle (p (see section 
wich was the zenithal emission angle of Cherenkov photons. 



3.5 I 
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Figure 14: Fitted shower asymmetry of electron induced showers for different pri- 
mary energies: Eq - lOGeV (top left) and Eq = 1 TeV (top right). The dots 
correspond to the mean expectation, if each shower is aligned before averaging as 
described in the text. The shown error bars represent the standard deviation of the 



individual bins. A parabola eq 16 was fit to the distributions. The black curve rep- 



resents the fit to the data and the red and green curves fits, where each bin content 
was increased or reduced by one standard deviation. The bottom left figure shows 
the range of the fitted uncertainties for different primary energy, and the bottom 



right figure the results of the parameters in eq 16 versus the primary energy. 



3. 7. Dependence of the angular distribution on the shower age 

The shape of the angular distribution of tracks dl/da and therefore the angular 
distribution of Cherenkov light has been found to be independent of the energy 
Eq. However, the electromagnetic cascade has an extension of a few meter (see 
section [33]). Within this evolution of the cascade it is plausible that large scattering 



angles occur more frequently at later stages of the shower than at earlier. Therefore 
we investigate how the angular distribution of tracks changes with the age of the 
shower. 



Figure 15 shows the angular distribution of the track length density versus the 
longitudinal length of the shower The distribution is found to be largely 
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zenith angle a [deg] 

Figure 15: Distribution of the relative track length versus the shower length z and 
the inclination angle or for a 1 TeV electron shower. The vertical color code corre- 
sponds to the histogrammed lengths I normalized per initial particle. 



dominated by the longitudinal evolution of the particle density and only a small 
difference of the angular distribution between the onset and the end of the cascade 
can be seen. 




zenith angle ct [deg] 



Figure 16: Angular distributions of the track length I for different slices of the 
longitudinal shower evolution. The figures show the histograms for a primary y 
with energy Eq = lOGeV (left) and Eq = I TeV (right) normalized to the track 
length per particle. The longitudinal distribution has been split into 3 slices of 
equal total track length. 



For a more detailed investigation we split cascades along the shower axis z 
into parts of different shower age. We chose three slices such that they contain the 
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same total track length and thus emit the same total amount of Cherenkov light.. 



The resulting angular distributions are shown in figure 16 Large differences are 



only seen in the very forward region a < 10°. With increasing shower age the track 
length becomes larger by about a factor 3. The differences for large angles a > 20° 
are comparably smaller, with about 10% change in yield and no obvious change in 
shape. The situation is similar for higher energy Eq. 




Figure 17: Angular distribution of the emitted Cherenkov photons for the different 
slices of the longitudinal shower evolution. The figures show the histograms for a 
primary y with energy £"0 - lOGeV (left) and Eq = \ TeV (right) normalized to 
the track length per particle. The longitudinal distribution has been split into three 
slices of equal track length. 



The resulting distributions of Cherenkov photons are shown in figure 17 The 
large opening angle of the Cherenkov cone leads to a substantial smearing of the 
angular distribution. Hence, the relatively strong effect into the forward direction 
for the angular track length distribution does not propagate to an equal strong vari- 
ation of the Cherenkov peak. Here, only an effect of +10% - 20% in the variation 
of the peak is observed. 

In summary, we find that the differential variation of Cherenkov light emis- 
sion along the length of the cascade is a relatively weak effect. A global angular 
distribution, as parameterized in section |3.5[ seems justified in particular when 
considering that the length of a cascade is short compared to the typical spacing of 
optical sensors in neutrino telescopes. However, for more detailed information we 



have repeated the angular parameterization eq 15 also for the three different slices 
in shower age separately. The results of the parameterization are given in table 



B.13 B.15 and B.15 in Appendix B 
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4. Summary and Conclusions 

We have simulated electromagnetic cascades with Geant4 for different primary 
particles and primary energies £"0. We have parameterized the total Cherenkov- 
light-radiating track length and its fluctuations, the longitudinal development of 
the cascade and the angular distribution of emitted Cherenkov photons. 

Our result for the total track length agrees within 2% with the result obtained in 
lfT2l but disagrees with other previous calculations. The relative size of fluctuations 
for different showers decreases oc 1 / V^o- 

The longitudinal profiles are found to be well described by a gamma distri- 
bution and a difference between e- and y is observed as expected [6|. However, 
quantitatively the position of the shower maximum deviates from the values in |f6l|, 
but agrees with EOl . For higher energies than lOOGeV we observe a change of 
slope in the elongation rate. 

The angular distribution of tracks in the cascade and the corresponding distri- 
bution of photons is found to be independent of £"0 and the type of primary particle. 

Systematic uncertainties of our parameterizations are related to the used spec- 
tral index and density of ice. They are of the order of 1 % for typical values of ice 
and can be corrected, e.g. by rescaling the observables of our parameterizations 
which depend on length scales such as the axis of the shower development to the 
correct media density and by using the correct Cherenkov angle. 

Also uncertainties of the used differential cross sections of electromagnetic 
processes in Geant4 can reach up to a few % but are generally substantially smaller 
ll2n . We do not expect these differential uncertainties to significantly affect our 
global results, particularly the total track length is not strongly affected. Hence we 
estimate a typical uncertainty of less than 1%. 

The LPM effect and dielectric suppression [22J are considered in the used ver- 
sion of Geant4 [23] but are only expected to become significant at larger energy 
than considered here. Electronuclear interactions have not been simulated, because 
their cross section is small. 

The fluctuations of the total track length for different individual showers are 
found to decrease with cx 1 / -^£^ and the related uncertainty is already smaller 
than 1 % for a few Ge V. 

For the determination of the uncertainties in the calculation of the angular dis- 
tribution of emitted Cherenkov light we investigate and parameterize the uncer- 
tainty related to azimuthal fluctuations and the evolution of the angular distribution 
with increasing shower age. The azimuthal asymmetry of tracks is found to be 
small (about 7% at 1 GeV) and to decrease oc 1/ V^o- When taking into account 
the large emission angle of Cherenkov photons this effect is expected to be further 
washed out for the angular distribution of Cherenkov light and is largely negligible 



19 



at high energies. Also the effect of the longitudinal shower evolution is small and 
results in differences of the width of the Cherenkov peak in the angular distribu- 
tion. This effect becomes even less important for distances larger than the scale 
length of the cascade and it is washed out e.g. by the scattering of photons when 
propagating through the medium. 

Appendix A. Geaiit4 configuration parameters used for tliis study 

In this chapter a summary of the defined media properties and physics pro- 
cesses is given. 



Appendix A.l. Materials 



Medium 


Density 


g 


Index of refraction 


Element 


Fraction of mass 


Ice 


0.91 


3 


1.33 


Hydrogen 
Oxygen 


88.81% 
11.19% 



Table A.l: Composition of ice as used in the Geant4- simulation. 



In Geant4 macroscopic properties of matter are described by G4Material and 
the atomic properties are described by G4Eleraent. A material can consist of mul- 
tiple elements and therefore represent a chemical compound, mixture as well as 
pure materials. For the performed simulations ice was used. Unless noted oth- 
erwise, the value n = 1.33 is used for the index of refraction. The simulated 



properties are summarized in table A. 1 



Appendix A.l. Physicslist 

All physics processes, which are used during the simulation must be registered 
in G4VUserPhysicsList. These simulations are based on the standard physics 
list G4EmStandardPhysics_option3. The included processes are summarized in 



table A. 2 The maximum energy for the cross section tables and the calculation of 



dEldx in Geant4 is lOTeV. 
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Particle 


Process 


Model 


r 


G4PhotoElectricEfFect 


G4PEEffectFluoModel 




G4ComptonScattering 


G4KleinNishinaModel 




G4GammaConversion 






G4eMultipleScattering 

G4eIonisation 
G4eBremsstrahlung 






G4eMultipleScattering 

G4eIomsation 

G4eBremsstrahlung 

(^Apnln c AviniViilati on 






G4MuMultipleScattering 

G4MuIomsation 

G4MuBremsstrahlung 

G4MuPairProduction 

G4MuNuclearInteraction 

G4CoulombScattering 




pi^,pi-,K*,K-,p* 


G4hMultipleScattering 
G4hIonisation 
G4hBremsstrahlung 
G4hPairProduction 






G4ionIonisation 
G4hMultipleScattering 
G4NuclearS topping 




all unstable particles 


G4Decay 





Table A. 2: Physics processes of most important particles used in the simulation. If 
no model is specified the default model is used. For hadrons and ions that are not 
listed multiple scattering and ionization are defined. 
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Appendix B. Parameterization results 





fit of /(£()) 




particle 


a/ cm GeV"^ 




e~ 


532.07078881 


1.00000211 




532.11320598 


0.99999254 


7 


532.08540905 


0.99999877 




fit of (TiiEo) 




e~ 


5.78170887 






5.73419669 


0.5 


r 


5.66586567 





Table B.3: Result of the parameterization of the effective track length versus pri- 
mary energy (section 3.2 1. The top table gives the results of fits of eqjTJfor the fits 
l(Eo) and the bottom the standard deviation CThatiiEo) of the fluctuations. 



particle 



a 



7 



2.01849 
2.00035 
2.83923 



1.45469 
1.45501 
1.34031 



0.63207 
0.63008 
0.64526 



Table B.4: Results of the energy dependence of the longitudinal fits, section 



3.3 



particle 


a/ sr ^ 


b 


c 


d/ sr-i 




4.27033 


-6.02527 


0.29887 


-0.00103 




4.27725 


-6.02430 


0.29856 


-0.00104 


y 


4.25716 


-6.02421 


0.29926 


-0.00101 



Table B.5: Averaged parameters describing the angular distribution of emitted 



Cherenkov light, section 3.5 
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particle energy / GeV a b 



e- 1 


1.96883 


0.62794 


3 


2.68228 


0.61705 


7 


3.30523 


0.64303 


10 


3.61481 


0.65247 


30 


4.16566 


0.62248 


70 


4.77945 


0.62823 


100 


4.98860 


0.63416 


300 


5.61779 


0.62033 


700 


6.10809 


0.62129 


1000 


6.24439 


0.61083 


3000 


7.03624 


0.63287 


7000 


7.60499 


0.64319 


10000 


7.86789 


0.65099 


e+ 1 


1.93375 


0.61535 


3 


2.74487 


0.63234 


7 


3.27811 


0.62541 


10 


3.55064 


0.64072 


30 


4.22803 


0.63283 


70 


4.77295 


0.62917 


100 


4.88106 


0.61512 


300 


5.55997 


0.61605 


700 


6.05207 


0.61644 


1000 


6.30800 


0.62317 


3000 


7.01259 


0.63265 


7000 


7.58980 


0.64096 


10000 


7.89891 


0.65069 


y 1 


2.49299 


0.60823 


3 


3.66575 


0.71860 


7 


3.99721 


0.66217 


10 


4.10746 


0.64095 


30 


5.08856 


0.68042 


70 


5.33660 


0.63890 


100 


5.60790 


0.64568 


300 


6.08445 


0.62277 


700 


6.61645 


0.63321 


1000 


6.78153 


0.62575 


3000 


7.47467 


0.64408 


7000 


7.96892 


0.64825 


10000 


23 8.13041 


0.65240 



Table B.7: Results of the fits of the longitudinal cascade development, section 



3.3 



particle 



energy / GeV 


a/ sr ' 


b 


c 


dl sr"^ 


3 


4.21013 


-6.01199 


0.30008 


-0.00099 


7 


4.21865 


-6.02752 


0.30088 


-0.00096 


10 


4.28687 


-6.02147 


0.29801 


-0.00108 


30 


4.29190 


-6.02932 


0.29836 


-0.00103 


70 


4.29202 


-6.02913 


0.29837 


-0.00104 


100 


4.24274 


-6.02275 


0.29962 


-0.00098 


300 


4.28351 


-6.02684 


0.29850 


-0.00104 


700 


4.28365 


-6.02672 


0.29850 


-0.00104 


1000 


4.28351 


-6.02679 


0.29849 


-0.00103 


3000 


4.28364 


-6.02685 


0.29850 


-0.00104 


7000 


4.28360 


-6.02690 


0.29850 


-0.00104 


10000 


4.28365 


-6.02689 


0.29850 


-0.00103 



e+ 3 


4.38344 


-t 


).03618 


0.29571 


-0.00113 


7 


4.19912 


-t 


) .00949 


0.30022 


-0.00095 


10 


4.24809 


-t 


) .01 802 


0.29919 


-0.00105 


30 


4.24807 


-t 


) .01 820 


0.29918 


-0.00104 


70 


4.28043 


-t 


).02557 


0.29851 


-0.00103 


100 


4.28069 


-t 


).02577 


0.29853 


-0.00104 


300 


4.28116 


-t 


) .02627 


0.29857 


-0.00105 


700 


4.28100 


-t 


) .02623 


0.29854 


-0.00103 


1000 


4.28106 


-t 


) .02620 


0.29854 


-0.00103 


3000 


4.28124 


-t 


) .02647 


0.29856 


-0.00104 


7000 


4.28125 


-t 


).02656 


0.29856 


-0.00103 


10000 


4.28134 


-t 


).02658 


0.29856 


-0.00103 


r 3 


4.15186 


-( 


) .02039 


0.30274 


-0.00085 


7 


4.14478 


-t 


) .00543 


0.30202 


-0.00093 


10 


4.30725 


-t 


).02368 


0.29745 


-0.00107 


30 


4.31395 


-t 


).03350 


0.29786 


-0.00103 


70 


4.25518 


-t 


).02358 


0.29922 


-0.00099 


100 


4.25584 


-t 


).02391 


0.29928 


-0.00101 


300 


4.25586 


-t 


).02383 


0.29926 


-0.00101 


700 


4.25605 


-t 


) .02399 


0.29926 


-0.00101 


1000 


4.28625 


-t 


).02812 


0.29849 


-0.00103 


3000 


4.28632 


-t 


).02803 


0.29849 


-0.00103 


7000 


4.28625 


-t 


).02805 


0.29848 


-0.00103 


10000 


4.28624 


-t 


).02800 


0.29848 


-0.00103 



Table B.9: Results of fits to angular distribution of emitted Cherenkov light, section 
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3.5 



particle energy/ GeV b (polmA cr 



e 


1 


0.07064716 


3.47699427 


0.02103594 




10 


0.02131804 


3.49345733 


0.00739830 




100 


0.00643688 


3.52309761 


0.00234943 




1000 


0.00202350 


3.48731664 


0.00074057 




10000 


0.00067586 


3.46473695 


0.00023550 


e+ 


1 


0.07112344 


3.48404544 


0.02099490 




10 


0.02097273 


3.51963176 


0.00724243 




100 


0.00663420 


3.48056236 


0.00232045 




1000 


0.00213402 


3.46339344 


0.00073926 




10000 


0.00065039 


3.50744179 


0.00023639 


7 


1 


0.06652718 


3.55998809 


0.02135593 




10 


0.02166958 


3.47688180 


0.00730116 




100 


0.00637762 


3.52556673 


0.00235480 




1000 


0.00207235 


3.48507985 


0.00074559 




10000 


0.00066903 


3.47379137 


0.00023858 



Table B.IO: Results of fits to the shower asymmetry, eq 16 



section 



3.6 



The used 



bin-size is 30° . The column (cr) gives the average standard deviation of fluctuations 
in each bin. 



particle 


/j/rad ^ 


e~ 


0.07029 




0.07064 


r 


0.06668 



Table B.ll: Results of the amplitude coefficient of the shower asymmetry, eq 



17 



section 3.6 



25 



particle energy / GeV a 



b 



c 



d 



e- 3 


142.9460887 


-9.57395706 


0.15172749 


-0.00266262 


7 


85.84316237 


-9.05308106 


0.16313946 


-0.00257812 


10 


77.91266880 


-8.91314499 


0.16406813 


-0.00284706 


30 


46.82551912 


-8.41395678 


0.17847633 


-0.00255785 


70 


38.02891569 


-8.19061093 


0.18434614 


-0.00257831 


100 


34.69834063 


-8.10546380 


0.18758322 


-0.00247138 


300 


28.97457908 


-7.91471330 


0.19335737 


-0.00247902 


700 


25.83417400 


-7.79936480 


0.19747241 


-0.00242657 


1000 


24.40584315 


-7.74365060 


0.19963483 


-0.00238802 


3000 


23.00296659 


-7.68029984 


0.20171506 


-0.00238785 


7000 


21.74138458 


-7.62373034 


0.20390902 


-0.00236181 


10000 


21.46673564 


-7.61084235 


0.20439010 


-0.00235443 


e+ 3 


137.9670566 


-9.54748394 


0.15276498 


-0.00256487 


7 


84.48637367 


-9.00950972 


0.16249928 


-0.00266916 


10 


71.63054893 


-8.84968017 


0.16696350 


-0.00265524 


30 


43.90262329 


-8.34648184 


0.18034774 


-0.00257627 


70 


39.84186034 


-8.23419693 


0.18270563 


-0.00258497 


100 


34.87615415 


-8.10856986 


0.18734870 


-0.00249897 


300 


29.07673188 


-7.91861470 


0.19324693 


-0.00248046 


700 


25.89164868 


-7.80387042 


0.19747634 


-0.00240558 


1000 


25.19834484 


-7.77379969 


0.19835006 


-0.00241144 


3000 


23.01752709 


-7.68236819 


0.20174971 


-0.00237442 


7000 


21.65536829 


-7.62068764 


0.20410342 


-0.00235278 


10000 


21.33114755 


-7.60423566 


0.20462732 


-0.00235435 



3 


69.70469246 


-8.90640987 


0.17075128 


-0.00211751 


7 


58.54718763 


-8.65787517 


0.17279750 


-0.00249754 


10 


60.52821188 


-8.64481073 


0.17011384 


-0.00280772 


30 


46.94065239 


-8.39738538 


0.17763706 


-0.00264893 


70 


33.61049111 


-8.06982820 


0.18848567 


-0.00247429 


100 


31.15779106 


-7.99239921 


0.19102320 


-0.00247773 


300 


26.69846316 


-7.83299913 


0.19630007 


-0.00244224 


700 


24.57774841 


-7.74984999 


0.19933205 


-0.00239464 


1000 


23.43128734 


-7.70044416 


0.20109158 


-0.00238882 


3000 


22.16347238 


-7.64277275 


0.20314619 


-0.00237112 


7000 


20.95984696 


-7.58781685 


0.20539395 


-0.00233737 


10000 


20.71724506 


-7.57557901 


0.20582216 


-0.00233654 
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to the angular distribution of Cherenkov light 



Table B.13: Results of fits of eq 
for the first third of the longitudinal cascade profile, section|3.7 
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particle 



energy / GeV 




a 




b 




c 




d 


a 

J 


1 

1 




J 


J 1 oyoDDZ. 


u 


^yjoyyjo^yj 


u 


UU IjD / ZU 


7 


2 


23159905 


-5 


66592786 





3870281 1 





00096165 


10 


2 


41790477 


-5 


69781676 





37321829 





00074460 


30 


2 


64370991 


-5 


72687457 





35830301 





00037598 


70 


2 


83551875 


-5 


77447099 





34890180 





00017227 


100 


2 


87505527 


-5 


78793324 





34732224 





00017138 


300 


3 


02408620 


-5 


81360770 





34020700 





00000261 


700 


3 


12021062 


-5 


83135934 





33600149 


-0 


00008677 


1000 


3 


14748514 


-5 


83561542 





33478931 


-0 


00011422 


3000 


3 


22395167 


-5 


85249577 





33188936 


-0 


00018386 


7000 


3 


28080118 


-5 


86421471 





32972990 


-0 


00022818 


10000 


3 


28101811 


-5 


86419504 





32973948 


-0 


00022808 



e+ 3 


2.04418647 


-5.57393706 


0.39850069 


0.00104086 


7 


2.19550814 


-5.63219580 


0.38766663 


0.00100450 


10 


2.38067758 


-5.68162469 


0.37507854 


0.00073911 


30 


2.71866151 


-5.75661338 


0.35520946 


0.00036425 


70 


2.74306888 


-5.75739258 


0.35363398 


0.00031489 


100 


2.85095217 


-5.77142273 


0.34768022 


0.00014361 


300 


3.01522146 


-5.81294569 


0.34066357 


0.00001562 


700 


3.10293302 


-5.82808660 


0.33673188 


-0.00006997 


1000 


3.14755994 


-5.83742082 


0.33491531 


-0.00011240 


3000 


3.21892957 


-5.85126920 


0.33206727 


-0.00017664 


7000 


3.26963124 


-5.86055633 


0.33005267 


-0.00022419 


10000 


3.30250950 


-5.86852133 


0.32892921 


-0.00024216 


y 3 


2.31187235 


-5.66208473 


0.37959321 


0.00077227 


7 


2.41282189 


-5.68451332 


0.37267906 


0.00068690 


10 


2.51707194 


-5.69019396 


0.36481000 


0.00052426 


30 


2.75219082 


-5.76825072 


0.35376221 


0.00034959 


70 


2.86679434 


-5.77857129 


0.34716252 


0.00015492 


100 


2.98152659 


-5.80586065 


0.34217497 


0.00003287 


300 


3.06538269 


-5.81769690 


0.33811151 


-0.00004781 


700 


3.16811393 


-5.84135741 


0.33410858 


-0.00013899 


1000 


3.19153423 


-5.84578346 


0.33314079 


-0.00015419 


3000 


3.26243729 


-5.86085963 


0.33046056 


-0.00020943 


7000 


3.30585081 


-5.86779646 


0.32870484 


-0.00025242 


10000 


3.30895722 


-5.87088305 


0.32877273 


-0.00024303 



Table B.15: Results of fits of eq 15 to the angular distribution of Cherenkov light 



for the middle third of the longitudinal cascade profile, section 3.7 
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particle 


energy / GeV 


a 


b 


c 


d 




J 




7. J / ZUJOJ'H- 




U.UUZDJoUl 




7 


85.99022343 


-9.05411822 


0.16307753 


-0.00257921 




10 


77.79962657 


-8.91179784 


0.16410414 


-0.00284619 




30 


46.88367741 


-8.41476975 


0.17842752 


-0.00255845 




70 


38.08621359 


-8.19192075 


0.18429858 


-0.00257872 




100 


34.72957609 


-8.10617168 


0.18754369 


-0.00247331 




300 


28.97315954 


-7.91468463 


0.19335970 


-0.00247882 




700 


25.95546681 


-7.80308554 


0.19724201 


-0.00243642 




1000 


24.48558230 


-7.74579464 


0.19946929 


-0.00239809 




3000 


23.04520886 


-7.68159582 


0.20162205 


-0.00239248 




7000 


21.71336645 


-7.62283497 


0.20397007 


-0.00235902 




10000 


21.46127003 


-7.61065501 


0.20440436 


-0.00235375 


e 


J 


1 ^8 1 nnsQ^fi 

1 JO. lUUoyju 


— y.j^o^ 1 '4^0 




—U.UUZ J U'4-oy 




7 


84.49824682 


-9.00981190 


0.16249033 


-0.00266739 




10 


71.69254483 


-8.85031418 


0.16693985 


-0.00265661 




30 


43.95573617 


-8.34754329 


0.18030217 


-0.00257905 




70 


39.85957512 


-8.23454455 


0.18268712 


-0.00258586 




100 


34.87521933 


-8.10854937 


0.18735001 


-0.00249890 




300 


29.07085848 


-7.91850474 


0.19325786 


-0.00247968 




700 


25.88875951 


-7.80377086 


0.19748135 


-0.00240549 




1000 


25.19809801 


-7.77379079 


0.19835066 


-0.00241142 




3000 


23.01750888 


-7.68236750 


0.20174975 


-0.00237442 




7000 


21.65536329 


-7.62068747 


0.20410343 


-0.00235278 




10000 


21.33114596 


-7.60423560 


0.20462733 


-0.00235435 


7 


J 


uy.oju / ju 


8 Qn8^^^Qn 
— o.yuo J J jyu 




A 009 1 1 Qlf^ 
— U.UUZl ly /u 




1 


58.28180667 


-8.65424155 


0.17295947 


-0.00248847 




10 


60.49942091 


-8.64435137 


0.17012934 


-0.00280746 




30 


47.11516858 


-8.40045586 


0.17750367 


-0.00265504 




70 


33.55526488 


-8.06827181 


0.18855456 


-0.00247353 




100 


31.16182031 


-7.99249181 


0.19101555 


-0.00247788 




300 


26.71881146 


-7.83341114 


0.19625992 


-0.00244509 




700 


24.56636625 


-7.74960945 


0.19935804 


-0.00239229 




1000 


23.43436187 


-7.70054100 


0.20108516 


-0.00238910 




3000 


22.16109501 


-7.64269786 


0.20315149 


-0.00237087 




7000 


20.96330326 


-7.58792493 


0.20538521 


-0.00233778 




10000 


20.71783561 


-7.57559703 


0.20582042 


-0.00233665 



Table B.17: Results of fits of eq 15 to the angular distribution of Cherenkov light 

for the last third of the longitudinal cascade profile, section 13. 7 
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